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Abstract 

In 1855 H. J. S. Smith [2j proved Fermat's Two Squares using the 
notion of pahndromic continuants. In his paper Smith constructed a 
proper representation as a sum of two squares of a prime number p, 
given a solution of + 1 = (mod p) , and vice versa. In this paper 
we extend Smith's approach to proper representations by sums of two 
squares in rings of polynomials on fields of characteristic different from 
2. Our approach will also work for other representations of integers, 
such as sums of four squares. 

We keep as far as possible the palindromic character of the rep- 
resentations. While our results are likely not new, we believe our 
extension of Smith's approach is new. 

Keywords: Fermat's two squares theorem; four squares theorem; continu- 
ant; integer representation. 
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1 Introduction 



Fermat's Two Squares Theorem, that is 

A prime number p is representable by the form + iff —1 is a quadratic 
residue modulo p. 

has always captivated the mathematical community. Equally captivating are 
the known proofs of such a theorem; see, for instance, [ll[2l[71[lTl[T3] . Among 
these proofs we were enchanted by Smith's elementary approach 0, which 
is well within the reach of undergraduates. We remark Smith's proof is very 
similar to Hermite's [7], Serret's [TT], and Brillhart's pQ. 

Two main ingredients of Smith's proof are the notion of continuant (Def- 
inition [2] for arbitrary rings) and the famous Euclidean algorithm (stated for 
arbitrary Euclidean rings in Table [T]) . 

Let us recall here, for convenience, a definition taken from [8| p. 148] 

Definition 1 Euclidean rings are rings R with no zero divisors which are 
endowed with a Euclidean function N from R to the nonnegative integers such 
that for all ri,T2 & R with ti ^ 0, there exists q,r G R such that T2 = qri + r 
and N(r) < N(ri). 

Among well-known examples, we are going to use the integers with N{u) = 
polynoms over a field with N(P) = 2'i^sree(P) ^^^^ ^(^q) ^ q_ 

The sequence (gi, g2, • • • , Qn) given by the Euclidean algorithm (Table [1]) 
on Ti and T2, with ti and T2 in R, is called the continuant representation of 
(ti, T2) as we have the equalities ti = [qi,q2, . . . ,qn]h (this notation is defined 
in next paragraph) and T2 = [q2, ■ ■ ■ ,qn]h unless T2 = 0. If T2 7^ 0, then h 
is a gcd of (ri,r2), else h = ti, in other words Rti + Rt2 = Rh, where Rt 
denotes the left ideal generated by r. 

Definition 2 (Continuants in arbitrary rings) Let Q be a sequence of 
elements (gi, q2, ■ ■ ■ , qn) of a ring R. We associate to Q an element [Q] of R 
via the following recurrence formula 

[] = 1' = 9i> bi> 92] = qiq2 + 1, and 
[qi, g2, • • • , Qn] = [gi, • • • , qn-i\qn + [gi, • • • , qn-2] «/n > 3. 

The value [Q] is called the continuant of the sequence Q. 



2 



Table 1: Euclidean division 



input: Two elements Ti,T2 

of a Euclidean ring R 

with its Euclidean function N. 
output: a gcd of rl, r2 followed by 

a sequence {qi, q2, ■ ■ ■ , qn), possibly empty, of elements of R 

1 

while T2 7^ 

find qn, t such that ri = g„r2 + 1 with N(t) < N(r2) 
n <(— + 1 

r2 ^ t 
endwhile 

return ri followed by (gi, ^2, • • • , (Jn) 



Let p be a prime number of the form Ak + 1. Smith's approach [2] relies 
on the existence of a palindromic sequence Q = (gi, . . . , g^, g^, . . . , gi ) of 
even length such that p = [Q]. He then derives a solution for z"^ + 1 = 
(mod p) with 2 < z < p/2, namely [g2, . . . , gs, gs, . . . , gi]. On the other hand, 
from this solution one can retrieve the palindromic sequence by applying the 
Euclidean algorithm to p and z, and then p = + y"^ where x = [gi, • • • , g^] 
and y = [gi, . . . ,gs-i]. 

Brillhart's optimisation [1] on Smith's construction took full advantage 
of the palindromic structure of the sequence (gi, . . . , g^-i, g^, gs, g^-i, . . . , gi) 
given by the Euclidean algorithm on p and z. He noted that the Euclidean 
algorithm gives the remainders 

ri = [qi+2, qs-i, qs, qs, qs-i, • • • , gi] (i = i, . . . , 2s - l), and 
r2s = 

so, in virtue of Smith's construction, rather than computing the whole se- 
quence we need to obtain 

X = r^-i = [gs,gs-i, • • • ,gi] 

y =rs = [qs-i,...,qi]. 

In this case, we have y < x < ^Jp, Brillhart's stopping criterium. 

In this paper we study proper representations x^ + (that is, with x, y 
coprime) in some Euclidean rings via continuants. In Section [2] we study 
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some properties of continuants in arbitrary rings. Section [3] is devoted to 
study proper representations + y'^ in some Euclidean rings. We examine 
later some representations xx + yy using rings with an anti-automorphism 
(Sections m and [5]), keeping the palindromic (or quasi-palindromic) nature, 
up to multiplication by units, of the continuant. While the results presented 
here are not likely new, we believe our presentation is new. 



2 Continuants 

In this section we derive some properties of continuants from Definition [21 
which we will refer to as Continuant Properties. 

P-1 The first property is the so-called "Euler's rule" [31 p. 72]: once we have 
all the products of subsequences of Q obtained by removing disjoint 
pairs of consecutive elements of Q, the continuant [Q] is given by the 
sum of all such products. The empty product is 1, as usual. 

P-2 If in a ring R we find a unit r commuting with all g^'s, then 



r -1 (-1)'= (-1)" 1 



[gi,...,g„] if n even 
r"^[gi, . . . ,g„] if n odd 



P-3 [gi, . . . , g„] = [gi, . . . , gi_i][gi+2, . . . , gn] + [gi, • • • , gi][gi+i, • • • , gn] 

To obtain this equality, it suffices to divide the products of subsequences 
oi Q = (gi, g2, • • • , gn) obtained by removing disjoint pairs of consec- 
utive elements of Q into two groups, depending on whether the pair 
g^gi+i {1 <i < n) has been removed or not. 

P-4 From the previous points follows 

[-9/1, -qh-u • • • , -gi, 0, gi, g2, . . . , g„] = 

[g.^2,g.+3,...,gn] for0</.<n-2 

1 li h = n — 1 

ii h = n 

P-5 [gi, . . . , g„] and [gi, . . . , gn-i] are coprime. Note that from the previous 
points P-[3] and P-[l| we have more precisely 

[-g„_i,-g„_2,-gi,0][gi,...,g„] + [-g„-i, -gn-2, -gijfe, • • • , gn] = 1- 
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2.1 Continuants in commutative rings 

If the ring R is commutative, then we have some additional properties. 

P-6 [gi,g2, •••,?«,] = [g„,. . . ,g2,gi]- 

P-7 The continuant [gi, . . . , g„] is the determinant of the tridiagonal n x n 
matrix A = (aij) with Oj j = qi for 1 < i < n, aj^j+i = 1 and Oj+i^j = — 1 
for 1 <i < n. 

The following identity due to Lewis Carroll (alias Charles Lutwidge Dodg- 
son) plays an important role in our study of continuants. 

Lemma 1 (Lewis Carroll Identity) Let C be an n x n matrix in a com- 
mutative ring. Let Cj^^...^j^;j^^...j^ denote the matrix obtained from C by omit- 
ting the rows ii, . . . ,is and the columns ji, . . . ,js- Then 

det(C) det(aj;i,i) = det{Ci,i) det{Cj.j) - det{Ci.j) det{Cj,i) 

where the determinant of the Ox matrix is 1 for convenience. 

The use of Lewis Carroll Identity and property P-[7] provides more prop- 
erties. 

P-8 [gi, g2, • • • , gnjfe, • • • , qn-i] = [gi, • • • , gn-i][g2, •••,?«] + (-1)" (when 
n > 2). 

P-9 In the case of even n with qi = g„+i_j for 1 < i < n, in other words if 
the sequence is palindromic, we see 

[gi, • ■ ■ ,gn/2,gn/2, ■ ■ ■,q2f + l = 

[Qi, • • • , C[n/2, qn/2, ■ ■ ■ , ^ijfe, • • • , qn/2, 9^/2, ■ ■ ■ , ^2] = 

iku • • • , (ln/2? + [gi, • • • , qn/2-i?){[q2, • • • , qn/2? + fe, • • • , qn/2~i?) 

More properties and proof techniques valid in the commutative case are 
given in [6l ch. 6.7] 
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2.2 Quasi-palindromic sequences 

Here again the rings are not necessarly commutative. 



Definition 3 An anti-automorphism of a ring R is an involution r ^ r 
such that T + a = T + a and Ta = ar for all elements t, a of R. 

Definition 4 Let R be a ring endowed with an anti- automorphism t ^ t. 
A quasi-pahndromic sequence of length n satisfies qi = qn+i-i for 1 < i < n; 
in particular, if n is odd the element q(^n+i)/2 satisfies q{n+i)/2 = Q'{n+i)/2- 

We have an obvious relation, 



P-10 ...,qi] = [qi, . . . g„] 

and counterparts of the properties P-[5] and P-O 

Lemma 2 (Noncommutative Lewis-Carroll-like Identity) Let t ^ t 

he an anti- automorphism in a ring R, satisfying moreover the conditions 



Let {qi, ■ ■ ■ qn) he a quasi-palindromic sequence of length n > 2 in R. The 
following relation holds 



Proof. We proceed by induction on n. Our basic cases are n = 2, 3. The 
result is clearly true for n = 2. 

For n = 3, since q2 is in the centre of R and qi commutes with q^, from 




if T = T then r helongs to the centre of R. 



(1) 



[Qi,---,qn][q2,---,qn-i] = [q2,---,qn][qi,---,qn-i] + {-IT 

= [gi,...,g„_i][g2, + (-I)''- 



[qi, Q2][q2, ga] - 1 = iqiq2 + I)(g2g3 + 1) - i 
= qiq2q2q3 + qiq2 + q2q?, 

= qiq2q2q3 + qiq2 + ^2^3 



we obtain qiq2q2q3 + qiq2 + q2q3 = fe, q3][qi, ^2] - 1 = [qi, q2, q3]q2- 
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For larger n, write E = [q2, . . . , qn-i] and F = [q^, . . . , g„_2]. Thus, E 
and F belong to the centre of R, and the following results come from the 
definition of continuant and Property P-[3] 

[gi,...,gn-i][g2, = {qiE+ [q3,...,qn-i]){Eqn + [g2, • • • , gn-2]) 

= qiE'^qn + qiE[q2, • • . , qn-2] + [Qs, qn-i\Eqn 
+ fe, • • • ,gn-i][g2, • • • ,gn-2] 

[gi, g2, • • • , gn-1, gn]^ = (gife, • • • , gn-i, gn] + [gs, • • • , g™])^ 

= {QiiEqn + [g2, • • • , gn-2]) + [gs, • • • , gn-ijgn + i^) ^ 

= qiEqnE + gi[g2, . . . , gn-2]-E + [gs, . . . , qn-i]qnE + FE. 

First note that [g2, . . . , g^] [g i, • • • ,gn-i] = fei, • • • , ^n-i] [g2, • • • , gn] because 
of the equality [gi, • • • , g^-i] = [g2, • • • , gn]- 

Since E commutes with the whole i?, we have 

qiE'^qn = qiEqnE 
qiE[q2, qn^2] = gi[g2, • • • , gn-2]^, and 
[gg, . . . , qn-i]Eqn = [gs, . . . , qn-i]qnE 

It only remains to check 

EE = [g2, . . . , gn-2] [gs, ■ ■ ■ , qn-i] + (-1)" 

= [^3, • • • , Qn-l] [g2, • • • , gn-2] + (-1)" 

but these equalities follows from the inductive hypothesis. □ 

Remark 1 For a quasi-palindromic sequence Q of length n > 3, we have 
[gi, g2, . . . , g„-i] = gi[g2, • • • , g„-i] + [gs, • • • , Qn-i] 

= gi[g2, •••,?«-!] + [g2,. . . ,gn-2] 



3 Proper representations in Euclidean rings 

As Smith's approach heavily depends on the existence of a Euclidean-like 
division algorithm, one may try to extend it to other Euclidean rings R. 
However, the uniqueness of the continuant representation may be lost. Ba- 
sically, the uniqueness of the continuant representation boils down to the 
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uniqueness of the quotients and the remainders in the division algorithm. 
This uniqueness is achieved only when i? is a field or R = W[T], polynomial 
algebra over a field F |9] (considering the degree as the Euclidean function). 
Note that in Z we guarantee uniqueness by requiring the remainder to be 
nonnegative. 

3.1 Non-commutative Euclidean rings 

We first use continuants to obtain a multiple zz + 1 of an element m of 
the form xx + yy, with x, y satisfying Rx + Ry = R and r i— t- r an anti- 
automorphism in the ring under consideration. 

Theorem 1 Let R be an Euclidean ring, and let t \-^t be an anti- automor- 
phism of R satisfying relations (QP- If m E R admits a proper representation 
m = XX + yy ( that is, with Rx + Ry = R), then the equation zz + lE Rm 
admits solutions. 

Furthermore, one of these solutions is equal to [gj, . . . , gT, gi, . . . , 
where (gi, g2, • • • , g^) is the sequence provided by the Euclidean algorithm on 
X and y. 

Proof. Let (x, y) (with N(x) > N(?/)) be a proper representation of m. 

If ?/ = then x is a unit, so m must be a unit and the ideal Rr is the whole 
ring R. Otherwise, the Euclidean algorithm on x and y gives a unit u and a 
sequence (gi, g2, • • • , g^) such that x = [gi, g2, . . . , qs]u and y = [g2, . . . , qs]u. 
Then 

XX = [gi, . . . , gs]nM[g7, . . . , gT], using Continuant Property P-fTOl 
xx = [gj, . . . , gT] [gi, . . . , qs]uu, since uu belongs to the centre of R 

yy= [gl, ...,g^][g2,...,gs]'uTZ 

m = XX + yy = [gj, . . . , gT, gi, . . . , qs]uu, using Continuant Property P-[3] 
Let z = [gj, . . . , g]", gi, . . . , gs_i], then applying Lemma |2] we obtain 

zz + l = (MM)"^m[g7II, . . . , gl, gi, . . . , g^-i] 

= (Mll)"^[gin, . . . ,gr,gi, . . .,qs-i\m 

since m is in the center of R 

That is, z satisfies zl + 1 E Rm, which completes the proof of the theo- 
rem. □ 
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3.2 Commutative rings: from x^^y^ to + 1 and back 

In this section we deal with the problem of going from a representation x^+y^ 
of an element m to a multiple + 1 of m and back. We begin with a very 
general remark valid in every commutative ring. 

Corollary 1 In a commutative ring R, if Rx + Ry = R then there exists 
some z & R such that + divides + 1. 

If R is Euclidean, we can explicit z and + with continuants. 

This relation can be interpreted using Lewis-Carroll Identity. The de- 
terminant of the tridiagonal matrix A associated to the palindromic se- 
quence . . . , gi, gi, . . . , g„) (see property P-[7] of continuants) is x'^ + 
with x= [gi, . . . , g„] and y = [q2, . . . , g„] if n > 1. 

Moreover (x^ + y'^){[qi, ■ ■ ■ g^-i]^ + [^2, • • • Q'n-i]^) = z^ + 1 where z is the 
determinant of matrix formed by the 2n — l first rows and columns of A (see 
properties P-[H]and P-fB]). 

A natural question is then: if m divides z"^ + 1 does there exist x, y such 
that m = + y^? But this question is much harder! We now give examples 
showing that no simple answer is to be expected. 

In general we cannot construct a representation x^ + y"^ of an element m 
from a solution of 2;^ -|- 1 = (mod m). 

As an illustration, consider the Euclidean domain F2[X] of polynomials 
on the field F2, where + 1 is a multiple of m = 2; + 1 for any polynomial 
z, square or not. Recall that in F2[X] the squares, and therefore the sums of 
squares, are exactly the even polynomials (i.e. the coefficient of X* is null if 
t is odd). Thus, the converse of Corollary [T] is false in F2[X]. 

Other examples are the ring Z[i] of Gaussian integers and its quotients 
by an even integer, since the squares and the sum of squares have an even 
imaginary part. Thus, no Gaussian integer with an odd imaginary part is a 
sum of squares, although it obviously divides = + 1. 

However, there are cases where the answer is positive 

Proposition 1 Let R be a commutative ring. If 2 is invertible and —1 is a 
square, say 1 + k'^ = 0, then x = (^) + (^) 

Proposition 2 Let R = ¥[X] be the ring of polynomials over a field F with 
characteristic different of 2 such that —1 is a non-square in F. 
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// m divides + with t coprime, then m is an associate of some 
+ with X, y coprime. 

Proof. We introduce the extension G of F by a square root u of —1. The 
ring G[X] is principal and z"^ + t^ factorises as {z — ut){z + ut). The two 
factors are coprime, since their sum and difference are respectively 2z and 
2ut, and 2 and u are units. Introduce gcd(m, z + ut) = x + uy, then x — uy 
is a gcd of m and z — ut owing to the natural automorphism of G. The 
polynomials x — uy and x + uy are coprime and both divide m. Thus, m 
is divisible by (x — ojy){x + uy) = + y^. On the other hand, m divides 
{z — ut){z + ut). Consequently, (x — uy){x + uy) is an associate of m. Since 
X — uy and x + uy are coprime, we have x, y are coprime. □ 

This gives some cases where a reciprocal of the first assertion in Corollary 
[H holds. 

Proposition 3 Let m be a non-unit of¥[X] and a divisor of z'^ + 1 for some 
z G F[X] with deg(2;) < deg(m). 

If¥ is afield of characteristic different from 2, where —1 is a non-square, 
then continuants provide a method for representing m as a sum of squares. 

Specifically, the Euclidean algorithm on m and z gives the unit u and 
the sequence {uqs,u~^qs-i, . . . ,u^~^^" qi,u^~^^'' qi, . . . ,u~^qs) such that x = 
[gi,...,g,] andy = [q2,...,qs]- 

Proof. Having a divisor m oi z^ + 1, we already now from Proposition [2] 
that the degree of m is even. We may assume that degree (2;) < degree (m) as 
we may divide z hj m. 

As m = (x^ + the Euclidean algorithm on some x and y will give the 
unit 1 and a sequence (gi, . . . , g^) such that x = [gi, . . . g^] and y = [g2, . . . , g^]. 
We may also assume degree(x) > degree(|/), otherwise, if x = \y -\- z with 
A a unit and z of degree smaller than the degree of x and y, then m = 

(((I + X^)y + Xz)' + 

■ ^ result we consider only the case where all 

gj's have degree at least 1. 

We then apply the Euclidean algorithm (Tabled]) to m and z, and obtain, 
by virtue of the uniqueness of division in polynomials, a sequence whose last 
non-null remainder is u. Consequently, m/u = x"^ + y"^ (see Property P-[2]of 
continuants). □ 

Below we illustrate this proposition through some examples. 
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First take m = 2X'^-2X^ + 3X^-2X + 1, then m divides {2X^ + Xy + 1. 
The Euchdean divisions give successively 



- 2X^ + - 2X + 1 ={2X^ + X)(X - 1) + - X + 1 



Here we have m/2 = [2 ■ (X - l)/2, 2'^ ■ (2X + 1), 2 ■ (2X + 1), 2-^ ■ (X - 1) /2] 
with M = 2, which gives m/2 = (X^ - X/2 + 1 /2f + (X/2 - 1/2)2 ^ ^2 ^ 
Since 2 is also a sum of two squares, we obtain m = {x + yY + {x — y)"^ — 
X^ + {X^ -X + lf. 

We find other examples among the cyclotomic polynomials. The cyclo- 
tomic polynomial $4„ e Q[X] divides X^^ + l. Thus is, up to a constant, 
a sum of two squares. Since $4n(0) = 1, the constant can be chosen equal to 
1. For an odd prime p, it is easy to check 



$4p(^) = j](-i)'=x2'= = j2 (-1)'=^"' + K (-i)'=x' 



For the small composite odd number 15, the computation gives 



$6o(^) = + - X^o - X« - X^ + X^ + 1 

= [-^7 -X^^ — -X^, — -X, — -X, -X, -X, — -X, — -X, — -X, -X, -Xj 

-X'''^ — — — x^ — x^ x^ x^ — x^ — x^ x^ — x^ x^ x^ 

X = [X, -X, -X, X^ - X, X, X] 

= X« - X^ + 1 
y = [-X, -X,X3-X,X,X] 

= X^ + X^ - X^ - X 
^eo{X)^x^ + y^ 



At this stage the following remark is important. 

Remcirk 2 If a polynomial with integer coefficients is the sum of squares of 
two polynomials with rational coefficients, it is also the sum of squares of two 
polynomials with integer coefficients 



2X3 ^ ^ =(2X2 _ X + 1)(X + 1/2) + {X/2 - 1/2) 
2X2 - X + 1 =(x/2 - 1/2)(4X + 2) + 2 
X/2 - 1/2 =2(X/4- 1/4). 



p-i 



k=0 
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For example, we see that SOX^ + UX + 1 = (5X + 3/5)2 ^ ^ 4/5^2^ 
it is also X2 + (7X + l)^. 

A proof was given by Gjergji Zaimi [12]. Another proof can be seen 
in m Sec. 5]. 

Remark 3 (Algorithmic considerations) To accelerate the computation 
in Proposition\^ we can resort to Brillhart's fJ]/ optimisation and stop when 
we first encounter a remainder r^-i with degree at most degree(m)/2. This 
will he the {s — l)-th remainder. In this context 

{Tg^i for odd s 

u fs-i for even s 

This observation follows from dividing m/u = [qs, . . . , gi, gi, . . . , g^] by z = 
[qs-i, . . . , gi, gi, . . . , gs] using continuant properties. 




[uqs, u ^qs-i, . . . , g2] for odd s 

u~^[uqs,u~^qs-i, . . . ,u^~^^" ^g2] for even s 



4 Four squares theorem 

We slightly generalise the formula for products of sums of four squares; see 
[TU| p. 135] (this was already known to Euler, see jS] p. 277]). 

Lemma 3 (Product formula) Let R be a commutative ring endowed with 
an anti- automorphism. Let x,y, z,u be elements of R. Then 

{xx + yy) {zz + uu) = {xz — yu) {xz — yu) + {xu + yl) {xu + yz) 

Proof. This can be seen by looking at the determinants in the equality 

xz — yu xu + yz 
—xu + y'z xz — yu 

□ 

The four squares product formula is the application of the lemma [3] to the 
case where R is the ring of Gaussian integers, with its conjugation. 

This product formula allows to reduce the proof of the four squares the- 
orem to the case of primes. 

We recall that in Z/pZ the element —1 is either a square or a sum of two 
squares; see [IHl p. 133] (this also was known to Euler [51, p. 279]). Thus, 
it suffices to prove that if a positive number m divides zz + 1 with z being 



X y 

-y X 



z u 
—u z 



X y 
-y X 
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a Gaussian integer, then m is also xx + yy with x, y both being Gaussian 
integers. 

By reducing z modulo m, we may assume \z\ < m/ a/2, and thus zz + 1 < 
(if m = 2, a parity argument shows the inequality remains valid). Here 
\z\ denotes the complex norm of z. 

Then, we have a succession of s equalities mirrii+i = ZiZl + 1 and Zi = 
qi+inii^i+Zi^i, where the sequence of positive integers m = mo, mi, . . . ,ms = 
1 is decreasing. At the end we have ms-ifrLs = z^-i^s-i + 1 and = z^-i- 

Now the quasi-palindromic sequence Q = (gi, g2, • • • , ?2,^), where the 
central pair is g^, gj if s is odd and gj, g^ if s is even, satisfies [Q] = m. Thus, 
we have a representation oim = xx + yy with x and y being the continuants 
of the sequences formed by the first s terms and first s — 1 terms of Q, 
respectively. The alternation of quotients and their conjugates is justified by 
Remark [TJ 

Consider the following example, where mo = 431 and zo = 54 + lOi. 

431-7 = (54 + 10i)(54-10i) + l ^ 54 + lOi = (8 + i)7 + (-2 + 3i) 
7-2 = (-2 + 3i)(-2-3i) + l ^ -2 + 3i = (-l + i)2 + i 
2-1 = (i)(-i) + l ^ i = i-l 

Hence Q = (gi, g2, gs) = (8 + i, -1 + i, i) and 

431 = [8 + i, -1 - 2, z, -z, -1 + z, 8 - z] 

= |[8 + i,-l-i,i]|2 + |[8 + i,-l-2]|2 = |17-5i|2 + I - 6 - g^l^ 
= 17^ + 52 + 6^ + 92 



5 Some forms representing integers 

Using the techniques of Section H] we may build other forms representing all 
positive integers. One such form is x^ — xy + y'^ + — zu + . 

Proposition 4 Each positive integer has the form x"^ —xy + y'^ + z^ —zu + v? 
with X, y, z, u integers. 

Proof. We note that v"^ — vw + w'^ is the norm of f + wj in the ring 
of Eisenstein integers where j = exp(227r/3). The same arguments as in 
Section H] reduce the task to primes and prove that each prime is either of 
the form zz or divides some zz + 1. 
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Now, if an integer m divides zl + 1, the division process provides an 
deterministic algorithm to find a representation m = xx + yy. Here again 
we reduce z modulo m and assume z'z < 3m^/4. Thus, we only have to be 
careful if rus-i = 2 to avoid the trap 2 ■ 2 = (1 — j)(l — j) + 1 by choosing a 
convenient quotient qs-i- □ 

Consider the following example where we try to represent mo = 40. We 
note 80 = 40 ■ 2 = (7 - 3j)(7 - 3j) + 1. With the quotient = 3 - j we 
would get 2-2 = (1— j)(l — j) + l. However, with the quotient gi = 3 — 2j, 
we get 2 ■ 1 = (1 + j)(l + j) + 1 and q2 = I + j- Hence 

40 = [3 - 2j, 1 + j][7 - 3j, 1 + j] + [3 - 2j][3 - 2j] 
= (-1 - 5j)(-l - 5j) + (3 - 2j)(3 - 2j) 
= 5^ - 5 • 1 + 1^ + 3^ + 3 ■ 2 + 2^ 
= 21 + 19. 



Corollary 2 Every positive integer has the form x^ + Zy^ + + 3m^. 

Proof. By Proposition H] we only need to prove that — + y^ has the 
form "i]? + g^. Indeed, 

• If X is even, say x = 2t, then x'^ — xy + y'^ = At^ — 2ty + y'^ = 3t^ + (y — t)^ 

• If y is even, say y = 2t, then x^ — xy + y"^ = 3t^ + (x — t)^ 

• If X and y are both odd, then x^ — xy + y"^ = ((x + ?/)/2)^ + 3((?/ — x)/2)^ 

□ 

Proposition 5 Each integer has the form x^ — Sy"^ + z"^ — 3m^. 

Proof. This follows reasoning as in Proposition HI The necessary ring is 
Z[-\/3] endowed with its natural anti-automorphism. □ 

In the following example we try to represent 19 and —19, noticing that 
19 . 2 = 72 - 3 ■ 22 + 1. 

19-2 = (7+y3)(7-y3) + l 

gi = 3 + 73 

2-1 = (1 + 0V3)(1-0V3) + 1 

ga = 1 + 0^3. 
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Hence 



19 = [3 + V3, 1 - 0V3][3 - V3, 1 + 0^3] + [3 + V3][3 - V3] 

= (4 + y3)(4 - 73) + (3 + V3){3 - V3) 
= 16-3 + 9-3. 

Then, to represent —19, we use -1 = 1 ■ 1 + (1 + y/3){l - V3) and the 
product formula (Lemma [3]) to get 

-19 = ((4 + y3)(l + VS) + (3 + V3)){A + y3)(l + \/3) + (3 + V3) 

+ ((4 + 73) - (3 + y3)(l - y3))(4 + 73) - (3 + y3)(l - Vs) 
= (10 + 6y3)(10 - 6^3) + (4 + 3y3)(4 - 3^3) = -8 - 11. 

Remark 4 We have proved the existence of decompositions of positive num- 
bers as sums of two norms of Gaussian integers or Eisenstein integers. The 
number of representations of positive numbers by these forms are given by Ja- 
cobi's theorem [?, Theorem 9.5] and Liouville's theorem [?, Theorem 17.3], 
respectively. 
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